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We use high-resolution inelastic neutron scattering to study the low-temperature magnetic ex- 
citations of electron-doped superconducting Pro.88LaCeo.i2Cu04-(5 {Tc = 21 it 1 K) over a wide 
energy range (4 meV< hu < 330 meV). The effect of electron-doping is to cause a wave vector 
(Q) broadening in the low-energy (huj < 80 meV) commensurate spin fluctuations at (7r,7r) and to 
suppress the intensity of spin- wave- like excitations at high energies (huj > 100 meV). This leads to 
a substantial redistribution in the spectrum of the local dynamical spin susceptibility x'i^)^ ^ind 
reveals a new energy scale similar to that of the lightly hole-doped YB2CU3O6.353 {Tc = 18 K). 



PACS numbers: 74.72.Jt, 61.12.Ld, 75.25.+Z 



High-transition-temperature (high- Tc) superconduc- 
tivity in copper oxides occurs when sufficient holes or 
electrons are doped into the CUO2 planes of their in- 
sulating antiferromagnetic (AF) parent compounds [l|. 
Given the close proximity of AF order and supercon- 
ductivity, it is important to understand the evolution 
of magnetic excitations in the AF ordered parent in- 
sulators upon chemical doping to produce metals and 
superconductors, as spin fluctuations may play a cru- 
cial role in the mechanism of superconductivity |^] . For 
the undoped parent compounds, where AF order gives a 
diffraction peak at wave vector Q = (7r,7r) or (0.5,0.5) 
(Fig. la), spin waves at energies (hoo) below 60 meV 
found by neutron scattering show commensurate exci- 
tations around (tt, tt) because of the large AF nearest 
neighbor exchange coupling (Ji > 100 meV, Figs. la,lc, 
and 3a) |3: Ij? I^l • Upon hole-doping to induce metallic- 
ity and superconductivity, the low-energy spin fluctua- 
tions of La2-x(Sr,Br)a^Cu04 (LSCO) form a quartet of 
incommensurate peaks at wave vectors away from (tt, tt) 
[a, lD, y, y| that may arise from the presence of static or 
dynamic spin stripes |lO(. For hole-doped YBa2Cu306+a? 
(YBCO) with X > 0.45, the magnetic excitation spectra 
have a commensurate resonance at (tt, tt) and incommen- 
surate spin fluctuations similar to that of LSCO below 
this resonance El E, E Q E 13113 . For energies 
above the resonance, the excitations are spin-wave-like 
for X < 0.5 I15LI17I I and become a "box-like" continuum 
at X = 0.6 ^^. In the extremely under doped regime 
{x = 0.353, Tc = 18 K), Stock et al [l^ showed that 
spin fluctuations are commensurate around (7r,7r) and 
have a damped spin resonance around 2 meV. 

While the evolution of spin excitations in hole-doped 
superconductors has become increasingly clear, it is cru- 
cial to determine the evolution of spin excitations in 
electron-doped materials, as particle-hole symmetry is 
an important ingredient of any theory purporting to ex- 



plain the mechanism of high- Tc superconductivity. Un- 
fortunately, due to the difficulty of growing large high- 
quality single crystals required for inelastic neutron scat- 
tering experiments, there exist only a few studies ex- 
ploring the low-energy {hoo < 10 meV) spin dynamics 
in electron-doped materials such as Ndi.85Ceo.i5Cu04 
(NCCO) [I2II2II: cind consequently the overall magnetic 
response of electron-doped materials remains largely un- 
known. Recently, we began to systematically investi- 
gate the evolution of AF order and magnetic excitations 
as Pro.88LaCeo.i2Cu04_(5 (PLCCO) is transformed from 
the as-grown AF insulator into an optimally electron- 
doped superconductor {Tc = 25 K) without static AF 
order through an annealing process with a minor oxy- 
gen content S modification [2l|, |2^ |2^ . We chose first to 
study underdoped PLCCO {Tc - 21 K, Tat - 40 K) for 
two reasons. First, this material is between the as-grown 
AF insulator and optimally doped PLCCO and therefore 
has a larger magnetic signal than that of the optimally 
doped PLCCO [IllIiiEii. Second, Pr3+ possesses 
a nonmagnetic singlet ground state in PLCCO, different 
from the Nd^+ magnetic ground state in NCCO 25]. 

In this Letter, we report the results of inelastic neutron 
scattering measurements that probe the low-temperature 
(T = 7 K) dynamic spin response of electron-doped 
PLCCO {Tc = 21 ± 1 K) for energies from 4 meV to 
330 meV. We determine Q and co dependence of the gen- 
eralized magnetic susceptibility x"{Qj^) m- We find 
that the effect of electron doping into the AF insulat- 
ing PLCCO is to cause a wave vector broadening in the 
low-energy commensurate magnetic excitations at (7r,7r), 
consistent with that of the NCCO 19, 20]. At high ener- 
gies {hco > 100 meV), the excitations are spin- wave- like 
rings, but with a dispersion steeper than that of the un- 
doped Pr2Cu04 4] and with a significant reduction in 
the spectral weight of the local dynamical spin suscep- 
tibility x''{^) m- ^ comparison of PLCCO and lightly 
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FIG. 1: Pictorial representation of the dispersions of the 
spin excitations in the a) insulating Pr2Cu04 and b) super- 
conducting PLCCO. b) Magnetic susceptibility measurements 
of Tc's in arbitrary units for the seven samples. The inset 
shows unit cell with exchange couplings Ji and J2. d-g) one- 
dimensional cuts through the spin excitations at hoo = 200=bl5 
meV along the [T, 1], [1,1], [1,0], [0, 1] directions. 



doped YBCO {x = 0.353) [13 reveals that the energy 
scale for x^\u;) in both materials is at ~2 meV, lower 
than the ~18 meV for the optimally doped LSCO p. 

We grew seven high quality (mosaicity < 1°) PLCCO 
single crystals (with a total mass of 20.5 grams) using the 
traveling solvent method in a mirror image furnace. To 
obtain superconductivity, we annealed the as-grown non- 
superconducting samples in a vacuum (P < 10~^ mbar) 
at T = 765 ± 1°C for four days. While both ends of the 
same cylindrical shaped crystal were found to have identi- 
cal Tc's, there are small (±1 K) differences in Tc's for sep- 
arately annealed samples. Figure Ic shows magnetic sus- 
ceptibility measurements on all seven crystals used in our 
neutron experiments. They have an average Tc = 21 ± 1 
K. For the experiment, we define the wave vector Q at 
(Qx^qy^Qz) as {H,K,L) = {q^a / 2ti , qya / 2i: , q^c / 27i) re- 
ciprocal lattice units (r.l.u) in the tetragonal unit cell 
of PLCCO (space group M/mmm, a = 3.98, and c = 
12.27 A). The seven PLCCO crystals were co-aligned to 
within 1° in the [H, H, L] zone using HB-1/HB-lA triple 
axis spectrometers at the High-Flux-Isotope reactor. Oak 
Ridge National Laboratory. Our inelastic neutron scat- 
tering experiments were performed on the MAPS time- 
of-flight spectrometer with the incident beam parallel to 
the c-axis (T-direction) of PLCCO at the ISIS facility 
|lq . Four different incident beam energies of Ei = 40, 
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FIG. 2: S{Q,uj) in the {H,K) plane at a) huj = 4zbl, b) 71.5± 
3.5, c) ISOzblO, d) 200±15, and e) 240=bl5 meV. The incident 
beam energy Ei = 40, 115, 200, and 400 meV data have 
counting times of 18, 60, 44, and 76 hours respectively with a 
source proton current of 170 /xA. f-j) Narrow Q-cuts passing 
through (0.5,0.5) along the dashed line directions shown in 
a-e). Upper open triangles in j) show a thicker cut at Uco = 
315 di 15 meV along the [1, 1] direction (integrated from -0.2 
to 0.2 along the [1, 1] direction). Solid lines are the calculated 
one-magnon cross sections from the linear spin wave fit to the 
data with Ji = 162 meV and J2 = 0. Horizontal bars below 
each cut show the instrumental resolution. 



115, 200, and 400 meV were used, and the scattering was 
normalized to absolute units using a vanadium standard. 
Figure 2 summarizes images of neutron scattering in- 
tensity 6'(Q, uj) centered about (tt, tt) at T = 7 K in units 
of mbarns/sr/meV/f.u. without any background sub- 
traction. At the lowest energy {hco = 4 ± 1 meV) probed 
(Fig. 2a), the scattering consists of a strong peak cen- 
tered at (tt, tt) with some phonon contamination evident 
at larger wave vectors. A constant-energy cut through 
the image reveals a commensurate peak on a fiat back- 
ground (Fig. If). The peak is significantly broader than 
the instrumental resolution (horizontal bar) and gives a 
correlation length of ^70 A. Upon increasing energy, the 
peak at (tt, tt) broadens in width (Fig. 2b) and weakens 
in intensity (Fig. 2g). With further increase in energy 
to hu = 145 ±15 meV, the scattering becomes a spin- 
wave-like ring (Figs. 2c and 2h). One-dimensional cuts 
through Fig. 2d at hco = 200 ±15 meV along four dif- 
ferent directions (Figs. Id-g) confirm that the scatter- 
ing is indeed isotropic and symmetric around (tt, tt) like 
spin waves. With increasing energy, the ring continues to 
disperse outward until magnetic scattering is no longer 
discernible atfia; = 315±15 meV (Fig. 2j). 



Figure 3a summarizes the dispersion of spin excitations 
determined from the cuts in Figs. 2f-j. The dashed boxes 
show the positions of crystahine electric field (CEF) ex- 
citations arising from the Pr^+ rare earth ions in the 
tetragonal structure of PLCCO 25]. Compared to in- 
tensities of Pr^+ CEF levels, Cu^+ spin fiuctuations in 
PLCCO are extremely weak and cannot be separated 
from the strong Pr^+ scattering at certain CEF energy 
positions. For reference, figure 3b shows an energy cut 
along Q = (0.5, 0.5, L) for the Ei = 115 meV data in 
which it is clear that CEF intensities at hco ~ 20, 85 
meV are significant relative even to the incoherent elas- 
tic scattering. 

To estimate the strength of the magnetic exchange 
coupling, we consider a two-dimensional AF Heisenberg 
Hamiltonian with the nearest (Ji) and the next near- 
est (J2) neighbor coupling (Fig. Ic). Since the zone 
boundary spin fiuctuations sensitive to J2 were unob- 
servable (Fig. 2j), we set J2 = and determined that 
Ji = 162 ± 13 meV renders the best fit to the data 
for hco > 100 meV. The corresponding calculated one- 
magnon cross sections are plotted as the solid lines in 
Figs. 2h-2j, and the resulting dispersion relation is shown 
as the solid line in Fig. 3a. At high energies {huj > 100 
me V) , the calculated spin wave dispersion coincides fairly 
well with the data, but the value of Ji is considerably 
larger than the hole- (La2Cu04, Ji = 104, J2 = —18 
meV) 5] and electron- (Pr2Cu04, Ji = 121 meV) |^ 
doped parent compounds (Fig. 3a). Assuming as-grown 
insulating PLCCO has similar AF exchange coupling as 
Pr2Cu04 4], our data suggest that the high energy spin 
fiuctuations in electron-doped PLCCO disperse faster 
than the spin waves of the insulating compound. There- 
fore they are unlikely to arise from the weak static AF 
order in the material 21]. 

Although high energy spin excitations in PLCCO are 
spin- wave- like, the observed scattering for hco < 80 meV 
is substantially broader than those predicted by the lin- 
ear spin- wave theory (Fig. 3a). A cut at ficj = 8± 1 meV 
through the (7r,7r) point along the [1,1] direction con- 
firms this point (Fig. 3c). For energies below 20 meV, 
the dispersion of spin fiuctuations has a nearest neigh- 
bor coupling of Ji = 29 it 2.5 meV (dash-dotted line in 
Fig. 3a). Therefore, the dispersion of PLCCO can be 
separated into two regimes. For energies {4 < huj < 80 
meV), the excitations are broad and weakly dispersive. 
For hcj > 100 meV, the fiuctuations are spin- wave-like 
with Ji larger than that of the insulating parent com- 
pound. 

In addition to determining the dispersion of spin ex- 
citations in PLCCO, the absolute spin susceptibility 
x"{Q^^) measurements in Fig. 2 also allow us to cal- 
culate the energy dependence of the local susceptibility 
X'\uj), defined as j x"{Q..^)d^Q/ j d^Q 0, Q- Fig- 
ure 4a shows how x"{'^) varies as a function of Tiix) for 
electron-doped superconducting PLCCO. Similar to hole- 
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FIG. 3: a) The dispersion of spin excitations in PLCCO. 
Points connected by solid lines denote only one continuous 
peak centered around (tt, tt) and represent the FWHM of 
a Gaussian fit with the instrument resolution deconvoluted. 
Solid boxes show energy levels of Pr^+ CEF scattering. Solid, 
dotted, dashed, dash-dotted lines show dispersions from lin- 
ear spin-wave fits with various exchange couplings, b) Log 
vs. linear energy- cut averaged from H — 0.45 to 0.55 along 
the [1, 1] direction and from H — —0.05 to 0.05 along the 
[T, 1] direction. CEF and phonon contamination energies are 
marked by arrows, c) Q-cut along the [1, 1] direction with 
/icj = 8 =b 1 meV with a Gaussian fit in a blue dashed line and 
the calculated one-magnon cross section in solid red with J 
= 121 meV. 



doped materials |a, 13, ISl, llfl O llS l2i , electron-doping 
suppresses the spectral weight of spin fiuctuations in 
PLCCO at high (> 50 meV) energies. For energies below 
50 meV, x" i!^) increases with decreasing energy and does 
not saturate dXfiuo = A^l meV, the lowest energy probed 
on MAPS. Assuming that the x"{^) in crystals of MAPS 
experiments (Tc = 21 ± 1 K) is similar to that of the pre- 
viously studied T^ = 21 K PLCCO sample IHIilli^, 
we can normalize the low-energy magnetic response of 
the Tc = 21 K sample obtained on SPINS spectrometer 
at NIST to that of the MAPS data. The outcome, shown 
as inset of Fig. 4, reveals a new energy scale of 2-3 meV 
for superconducting PLCCO. 

We are now in a position to compare the spin exci- 
tations of electron-doped PLCCO with that of the hole- 
doped LSCO ^^ and YBCO [ll El ElM El • For 
hole-doped materials such as LSCO |aE,llJaland YBCO 
with X > 0.45 Q, EIIEES^ ^^^ low-energy spin fiuc- 
tuations are incommensurate and display an inward dis- 
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FIG. 4: Energy dependence of local susceptibility x\^) iii 
PLCCO determined from integration over wave vector of 
the observed magnetic scattering around (tt^tt) 0, |2g . The 
dashed line shows X {uj) x 1/5 for La2Cu04 |3- Since Q-cuts 
were made along the [1, 1] direction, the background scatter- 
ing can be approximated by a constant. The blue triangles 
indicate normalized cold neutron triple- axis data on a Tc = 21 
K PLCCO at energies below 3 meV obtained on the SPINS 
spectrometer at NIST center for neutron research 23]. The 
inset shows an expanded view of x\^^ vs liuj at low ener- 
gies. The solid line is a damped Lorenztian on a constant 
background. 



persion toward a resonance point with increasing energy. 
This is not observed in electron-doped materials. In- 
stead, spin fluctuations in PLCCO have a broad com- 
mensurate peak centered at (tt, tt) at low-energies (< 50 
meV) which disperses outward into a continuous spin- 
wave, ring-like scattering at high energies (> 100 meV), 
similar to lightly doped YBCO with x = 0.353 [ig. At 
present, it is not clear how theoretical models based on 
spin stripes [Ifl can reconcile the differences in spin exci- 
tations between the hole- and electron-doped materials. 

For hole-doped LSCO with Tc = 38.5 K, the mean- 
squared fluctuating moment (m^) = J x'\^)duj inte- 
grated up to 40 meV is (m^) = 0.062 ± 0.005 fi%f.u.-^ 
|9|. For comparison, (m^) calculated from x'\^) up to 
40 meV is only 0.024 ±0.003 /i|f.u.-^ in PLCCO, about 
three times smaller than that of LSCO. The total fluc- 
tuating moment integrated from to 300 meV (Fig. 4) 
gives (m^) = 0.089 ± 0.009 /i^f.u."^, a value an order of 
magnitude larger than the static moment squared (0.0016 
/i^f.u."^) [2l|. Since the total moment sum rule for spin- 
I Heisenberg model requires one-magnon fluctuating mo- 
ment squared to be smaller than the ordered static mo- 
ment squared [2^, the observed high-energy spin- wave- 
like excitations are unlikely to arise from the small or- 



dered moment. 

In the standard Hubbard model and its strong- 
coupling limit, the t-J model with only the nearest- 
neighbor hopping t, there should be complete particle- 
hole symmetry and therefore the electron- and hole- 
doped copper oxides should behave identically. The 
observed large difference between incommensurate and 
commensurate spin fluctuations in hole- jy, [3, ]^ Isj and 
electron-doped materials 19, 20, 21, 22, |2^ has mostly 
been attributed to their differences in the strength of sec- 
ond nearest-neighbor {t') and third (t'') nearest-neighbor 
hopping. This also explains their differences in Fermi 
surface topology within the t-J model J23, 1221 l^fll^ ^^~ 
though it may also be due to their proximity to two dif- 
ferent quantum critical points |2l|. In the most recent 
calculation using the slave-boson mean-field theory and 
random phase approximation |30^, incommensurate spin 
fluctuations at (0.37r, O.Ttt) have been predicted for opti- 
mally doped NCCO. However, this is not observed in our 
PLCCO (Figs. 2a-2e). Similarly, the energy dependence 
of the x^\Q,uj) at Q = (tt, tt) has been predicted to ex- 
hibit a peak between OAu/J 0, l^^ and OAu/J [2^. 
While qualitatively similar to the predictions, x''{^) ^^ 
Fig. 4 has a peak at a much smaller energy of 0.02a;/ J. 
Comparison of future calculations in absolute units with 
our data should determine whether itinerant magnetism 
models can account quantitatively for the observed dy- 
namic susceptibility in superconducting PLCCO. 
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